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Abstract 

We study partial differential equations of second order (in time) that possess a hierarchy of infinitely 
many higher symmetries. The famous Boussinesq equation is a member of this class after the extension of the 
differential polynomial ring. We develop the perturbative symmetry approach in symbolic representation. 
Applying it, we classify the integrable equations of 4th and 6th order (in the space derivative) equations, as 
well as we have found three new 10th order integrable equations. To prove the integrability we provide the 
corresponding bi-Hamiltonian structures and recursion operators. 



1 Introduction 

The variety of integrable evolutionary equations, i.e. equations explicitly resolved with respect to the first time 
derivative 

^\^3 ^X , I^XX] '"3 ^XX...x) ) 

have been extensively studied and most comprehensive results of classification of integrable equations have been 
obtained for this class of equations (see for example PE])- 

In this paper we study necessary conditions for the integrability of non-evolutionary equations of the form 

u tt = adlu + f3d%u t + K{u,u x ,u xx , ...,dl~ l u,u u u tx ,Utxx, --^dl^Ut) , p > q , a,/3eC (1) 

and analyse equations satisfying these conditions. We call an equation to be integrable if it does possess an 
infinite hierarchy of higher symmetries. The famous Boussinesq equation 

u tt — u xxxx { u )xx (2) 

belongs to this class (after the extension of the differential polynomial ring). Recently all integrable equations 
of the form 

<M;t — l^xxx 

+ F(u,u x ,u xx ,ut,u t x) (3) 
have been classified and comprehensively studied in 0]. 

Equations of the sixth order (p = 6) of the form JQ) have been studied in connection with reductions of the 
Sato hierarchies corresponding to KP, BKP and CKP equations [SJ|S]. Sixth order equations and their Lax 
representations can also be derived from the systems studied by Drinfeld and Sokolov [7||S]. 

The aim of our paper is to look in depth in the structure of hidden higher symmetries, to develop a method that 
would allow to solve the classification problem for integrable equations of relatively low order. It is a natural 
development of the Perturbative Symmetry approach |H] based on the standard symmetry approach and 
symbolic method developed in |l(Jj . In the spirit of perturbation theory we define approximate infinitesimal 
symmetries and approximate recursion operators. Equations that have approximate recursion operator possess 
an infinite hierarchy of approximate symmetries and it is natural to call such equations approximately integrable. 
A similar idea of approximate integrability, based on approximate symmetries, conservation laws and recursion 
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operators was discussed in the context of obstacles to asymptotic integrability |1 11112] . In this paper we formulate 
explicit and easily verifiable criteria for approximate integrability. The cases of odd and even order p equations 
are essentially different. The general theory presented in Section 2 is suitable for the both cases, but in this 
paper we apply the theory to even order equations. The results of the odd order equations will be published 
somewhere else. 

In Section 3 we show that the method proposed in this paper does reproduce all known integrable cases of 
equation (|TJ (for p = 4, 6) as well as enable us to find three new 10th order equations. We present recursion, 
Hamiltonian, symplectic operators and the corresponding Lax representations for integrable equations of 6th 
and 10th order. 



2 Approximate symmetries and recursion operators 

A non-evolutionary equation 

u tt = K(u,u x ,u xx , ...,d2u,u t ,utx,u t xx, —,d™u t ) , (4) 

can always be replaced by a system of two evolutionary equations 

u t = v, _ 
v t = K{u, u x , u xx , d%u, v, v x , v xx , 



Non-evolution equation may have other representations in the form of an evolutionary system of equations. 
If K — D X {G), where D x is a total x-derivative, then the system of evolutionary equations 

u t =v x , v t =G. (6) 

also represents @. 

For example, eliminating variable v from the following equations 
Case 1: 

u t = v , v t = u xxxx + {u 2 ) xx , (7) 



Case 2: 

u t = v x , v t = u xxx + {u 2 ) x , (8) 



Case 3: 

u t = v xx , v t = u xx + u 2 , (9) 

we receive the same Boussinesq equation (J2J on variable u (variable v has different sense and scaling dimension 
for these three cases). 

There is a subtle, but important difference between the above representations of the Boussinesq equation J3J) • It 
is easy to verify that the system Q has only a finite number of local infinitesimal symmetries (i.e. symmetries 
whose generators can be expressed in terms of u, v and a finite number of their derivatives) while equations 
(|HJ) and have infinite hierarchies of local symmetries. The reason is simple - infinitesimal symmetries of 
equation JHJl, which depend on variable v, cannot be expressed in terms of u and its derivatives, since formally 
v = D^ut- There is only a finite number of symmetries of JHJ that do not depend explicitly on the variable v , 
and only these symmetries do survive the change of variables to equation (7J . In the extension of the differential 
ring by the element v = D~ x ut system (J7J) has infinitely many infinitesimal symmetries. 

These examples shows that one has to be rather careful with the definition of symmetries and more rigorous 
language is indeed required. 



2.1 The ring 1Z of differential polynomials. 

In what follows we assume that all functions, such as K and G in formula JSJ and ©, are differential polynomials 
of variables u, v, u x , v x , u xx , .... Instead of functions u, v and their partial x-derivatives we introduce an infinite 
sequence of dynamical variables {uq, Uq, u±, vi, 1*2, ^2, ■ • ■ } assuming the identification 

u =u, v = v, u n = d x u, v n = d x v. (10) 
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Often we will omit the zero index and write u and v instead of uq and vq. 

We denote TZ the ring of polynomials over C of infinite number of dynamical variables. We assume that 1 ^ TZ. 
Elements of the ring TZ are finite sums of monomials with complex coefficients and therefore each element 
depends on a finite number of the dynamical variables. The degree of a monomial is defined as a total power, 
i.e. the sum of all powers of dynamical variables that contribute to the monomial. Monomials are eigenfunctions 
of the following operator 

D = X U + X V (11) 

where 

fe>o OUk k>o OVk 

The eigenvalue of Dq is the degree of a monomial. We say that a polynomial is homogeneous of degree n if 
every its monomial is of degree n. The ring TZ has a natural gradation 

n = nn > nn ■ nm c nn+m > nn = if e 71 1 D of = n f) ■ 

nez+ 



Elements of 1Z 1 are linear functions of the dynamical variables, 7Z 2 quadratic, etc. It is convenient to define a 

^k>n 



'little-oh" order symbol o(TZ n ). We say that / = o(7Z n ) if / € 72. fc , i- e - the degree of every monomial of 



/ is bigger than n. 

Let v be two positive integers 1 , which we call the weights of u and v respectively and denote W(u) = 
fj,, W(v) — v. We say that a polynomial / 6 1Z is a homogeneous polynomial of weight A (and write W(f) — A) 
if X^f = A/, where 

X^ = (/xX„ + i/X t ,+X) , X= ^ k(u k ^-+v k ^-) . (12) 



fcez 



+ 



Eigenvalues of the operator X^ v are positive integers and the minimal eigenvalue is equal to min(/x, v). Mono- 
mials are obviously eigenfunctions of X^ v . The weighted gradation of 7Z is defined as 

TZ= K n , Tin ■ K m c K n+m , K n = {feTZ\ X„ v f = nf} . 

n£Z+ 

For example if fi = 1 , v = 2 then 

7Zi=C-u , 72.2 = spaxi(u ,ui,vo), TZs = span(u , uoUi,uqVq,U2,v-l), . . . , 

Operators Dq and X^ v commute and therefore we can define a degree- weighted gradation of the ring 1Z 

iz= iz;, hz-k™ck;+™, tz; = {/ en\x^f = P f , D f = nf}. 

Linear subspaces TZ n are infinite dimensional, subspaces 7Z n and 72.™ are finite dimensional (if we admit zero or 
negative weights, then the dimension of 7Z n would become infinite). 



The ring TZ is a differential ring with derivation 



^ = ^( Ufc+1 — + Wfe+1 — ) . (13) 



fc>0 

Derivation D x represents the x-derivative (c.f. 110(1 '). Since 1 ^ TZ, the kernel of the linear map D x : TZ 
ImD x C TZ is empty and therefore D~ l is defined uniquely on ImD^. It is easy to verify that 

[D ,D X ] =0, [X„ V ,D X ] = D X 

implying 

n • TP™ | > 1?" 



1 In general, the weights could be rational numbers including zero, but in this paper we consider only positive integer weights. 
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With any evolutionary system of equations 



u t = P(u,v,u x ,v x , . . .) , v t = Q(u, v,u x ,v x , . . .) , (14) 
where P, Q are differential polynomials, we associate an infinite dimensional dynamical system 

4j:(uk) =D%P(u ,vo,ui,vi, ...), ■^(yi t )=D*Q(uo,v ,ui,vi,...), k = 0,1,2,.... (15) 
It is sufficient to define the first two equations (u = uq, v — vq) of l|15|l 

j t ^ = P, j t v = Q> P,Qen, (i6) 

the rest equations of this infinite system can be obtained by the application of D x , k = 1, 2,3, . . . assuming 
that A4 = |A*. 

The system ()15(l provides a ^-derivation of the ring 1Z 



fc>0 



The derivations D x and D t commute, i.e. [D x ,D t ] = 0. Derivations of the ring 1Z that commute with D x 
we call evolutionary derivations. There is one to one correspondence between evolutionary systems of PDEs 
and evolutionary derivations. We say that evolutionary derivation D t is generated by (P,Q), therefore it is 
convenient to adopt notation D a where a = (P, Q)^ 1 denotes a vector column of the right hand side of (|1 6|) . 

The evolutionary system Ijl4(l and corresponding derivation D t are called homogeneous of weight s and write 
W(D t ) — s if P,Q are homogeneous elements of 1Z and W(P) — /z = W(Q) — v = s. If W(D t ) = s then 



Dt '■ TZn 1 ► TZn 



+s ■ 



Moreover, if P £ lZ k s+u , Q € K* +v then D t :Tl™i — > n ™+k-i 



For example Ao, cf. (jl 1|> and A, cf. (|13|l are a homogeneous derivations of 1Z generated by (uo,vo) and 
(wi,fi)^ r respectively with W(Dq) = 0, W(A) = 1, while X (|12|l is a non-evolutionary derivation of 72.. 

It follows from l|17|l that for any element g G 1Z we have 



/c>0 

where 



fc>0 Ai>0 

are two linear differential operators. The sums in l|18f) are finite, since g depends on a finite number of dynamical 
variables. Thus, with any g e 1Z we associate two linear operators, a vector row 

.9 -» .9* = {g*u , g*v) , (19) 

which is called the Frechet derivative of g. 

Let us consider a two dimensional liner space JC — 1Z over C of vector columns whose entries are elements of 1Z. 
Every element a = (P, Q) tr € JC generates an evolutionary derivation A (A cf. formula (TTJl) of the ring and 
vice versa, with any evolutionary derivation we associate a generator, i.e. an element of JC. Derivations of 1Z 
form a Lie algebra over C with the Lie product of two derivations defined as a commutator of the derivations. 
The evolutionary derivations form a subalgebra of the Lie algebra of all derivations. Indeed, it is easy to verify 
that if D a and are two evolutionary derivations generated by a and b respectively, then the commutator 
A o D\y — Z?b ° A is an evolutionary derivation A a ,t>] generated by the Lie bracket [a, b] s JC of elements a 
and b where 

[a,b] =b*(a)-a*(b), (20) 
or in components, assuming a=^ P ^,b=^^^we have 

r ii I S* U (P) + S* V {Q) — P*u(S) — P* V (T) \ (OT\ 

[a ' J ~ I T m {P) + T* V {Q) - Q* U (S) - Q* V (T) J ■ 
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In particular [D a ,Db] = [a, b] = 0. Wc shall assume that the linear space K. is equipped with the Lie 

bracket (|20|l and therefore it is a Lie algebra. 

The components of the Frechet derivative are examples of differential operators with coefficients from the ring 
1Z extended by the base field C: 

TZ C = C © K. 

We define formal series of the form 

A = ^2 a N-kD^~ k , a m e1Z c , togZ. (22) 

The order of the formal scries (|22|l is N (we assume that the leading coefficient a at ^ 0). The positive part 
of A, denoted by A + , which is a finite sum of the form A + = X^fcLo a N-kD x ~ k is a differential operator with 
obvious action on 7Z, while a formal series is not an operator, since no action on the ring JZ is defined. The sum 
of formal series is defined obviously, a multiplication (composition) is defined by 

a n DZ o b m D™ = E ( fc ) <^ D l iK)D7 +n - k ■ (23) 

For positive n the sum ((23(1 is finite since the binomial coefficients 

n \ _ n(n - l)(n - 2) ■ ■ ■ {n - k + 1) 
k ) ~ k\ 

vanish for k > n, for negative n the composition is well defined in the sense of formal series. Formal series form 
an associative ring and we denote it by lZc[D x ]. 

Suppose A is a formal series (differential operator) of the form 122(1 . then the conjugated formal series (differential 
operator) is defined as 

A^=J2(-l) N - k D?~ k oa N _ k . (24) 

We also consider formal series and differential operators with matrix coefficients, assuming that the entries of 
the matrix belong to IZc (the Frechet derivative 1(19(1 is an example of the operator with vector coefficients) . In 
this case the coefficients of the conjugated series 1(24(1 have to be transposed. 

For any g G 7Z we define a vector column of variational derivatives 5g 

*-(5S)-*!«-(tw) (25) 

There is a useful statement ^3], that 5g = if and only if g G lmD x , i.e. there exists h <E 1Z such that 
g = D x (h). Thus, the variational derivative is well defined on the factor space TZ/YvnD x , indeed two elements 
a,b G 1Z have the same variational derivative iff their difference is a total derivative, that is, a — b G Im_D x . 



2.2 Symmetries. 

There are many equivalent definitions of infinitesimal symmetries of PDEs (based on infinitesimal transforma- 
tions of solutions, commuting vector fields, etc., see for example PIE]). In this paper we adopt a definition, 
which is quite convenient for practical computations as well as for theoretical considerations. 

Suppose we have an evolutionary system ((16(1 and the column vector of the right-hand side is a = (P, Q)^ T . We 
say that b G K, is a symmetry (a generator of an infinitesimal symmetry) if [a, b] = 0. The evolutionary system 
corresponding to b is compatible with ((16(1 and often it is called symmetry as well. Due to the Jacobi identity a 
commutator of two symmetries is again a symmetry, therefore symmetries form a subalgebra of the Lie algebra 
K. 

Let us consider equation J2J which is represented by evolutionary system |[SJ. In this case a = {v,K)^ 1 . 
Assuming b = (S,F)^ T , and evaluating [a, b] = using 1(21(1 . we receive 

S*u(v ) + S. V (K) =F, F* u (v ) + F* V {K) = K* U (S) + K* V (F) , 

which can be rewritten as 

D t (S)=F, D t (F) = D b (K) , 
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where 

a = E + ^w^) . a, = E (^(^ + ^(^) (26) 

A symmetry generator b is completely determined by its first component 5, which satisfies equation 

D 2 t (S) = D b (K) . (27) 

Equation 1)27(1 together with 1)26)1 can be taken as a definition of symmetry for evolutionary system . 

Evolutionary system © represents equation (@J in the standard set of dynamical variables 1)10(1 if function 
K = D X (G) does not depend on the variable v (it may of course depend on v x — u t ,v xx — u tx , . . .). Indeed, the 
variable v cannot be represented in terms of u and its derivatives. Its symmetries may explicitly depend on v and 
therefore the corresponding infinitesimal transformations are non-local, i.e. cannot be expressed in terms of u 
and its derivatives. System (JHJ has infinitely many symmetries, but only a finite number of them does not depend 
on v explicitly. Thus, only a finite number of symmetries of the Boussinesq equation are local. The meaning of 
the notion "locality" depends on the choice of the set of dynamical variables. Representing the Boussinesq in 
the form JfjJ) we introduce a variable v and extend the standard set. In this extended set of dynamical variables 
the Boussinesq equation has infinitely many local infinitesimal symmetries. In the case of the system (0, 
there is another natural extension of the standard set of dynamical variables 1)10)1 . Due to the structure of the 
first equation ut — v\ we can introduce a "potential" w such that u = wi, u% = u>2, ■ ■ ■ Uk = Wfe+i- In this new 
extended set symmetries that depend explicitly on the variable v become local (v — w t ). Thus, the potential 
version of the Boussinesq equation 

w u — Wi + 2wiW2 (28) 

has infinitely many local symmetries. In the case of the Boussinesq equation one can extend the standard set 
one step further, introducing the next potential z such that z\ = w, Z2 — w\ — u, . . . z n+ 2 — uin+i = u n and 
therefore z u — z^ + (2:2) 2 • In the case of homogeneous equations introduction of a potential reduces the weight of 
the variable by one (in the case of the Boussinesq equation we have W(u) = 2, for the potential Boussinesq 
equation (3]jj) W(w) = 1, and W{z) = 0). 

For system © we have a = (vi, G)^ r . Denoting b = (5, F)* r it is easy to show that condition [a, b] = can be 
written as 

D t (S) = D x (F), D t (F) — D h (G) , (29) 



where 



a = e(^+^)^) < ^E^^+iT 1 ^*)^) . 



It follows from the first equation of the system 1)29)1 that S is either a density of a local conservation law of the 
system @ ora total derivative (S £ ImD x ). In both cases F = D~ 1 D t (S) £ 1Z and F can be found uniquely if 
S is known. Eliminating F from ()29)l we arrive to equation D^(S) — D\,[K), K — D X (G), which looks exactly 
as (127)1 . but the difference is in the definition of derivations 1)30)1 . 

Definition 1. We say that differential polynomial S GlZ is a generator of approximate symmetry of degree n 
of system 0) if 

D 2 t S - D h K = o{TZ n ) 

(D t and D\> are given by \3U\) ). We say that system has an infinite hierarchy of approximate symmetries of 
degree n if every member of the hierarchy is an approximate symmetry of degree n. 



Approximate symmetries of degree 1 are simply the symmetries of the linear part of the system. Any system 
(|6*|) has infinitely many approximate symmetries of degree 1. The requirement of the existence of approximate 
symmetries of degree 2 is very restrictive and highly non-trivial. An equation may have infinitely many ap- 
proximate symmetries of degree 2, but fail to have approximate symmetries of degree 3. Integrable equations 
have infinite hierarchies of approximate symmetries of any degree. The degree of the hierarchy of approximate 
symmetries is, in a certain sense, a measure of integrability. 

A homogeneous equation has homogeneous symmetries. Indeed, suppose S generates a symmetry and is a sum 
of components S = Si + ■ ■ ■ +Sk of different weights W(S P ) — n p , then it is easy to check that each homogeneous 
component S p generates a symmetry. We define the order of a homogeneous symmetry S as the weight of the 
corresponding evolutionary derivation ord(S') = W{D\ 3 ) = W(S) — W(u), which is obviously homogeneous. 
For example ord(ui) = W(D X ) = 1, for systems of equations 0, (JSJl and (0 we have ord(u t ) = W(D t ) = 2. 
Therefore the order of a homogeneous non-evolutionary equation of the form Q is natural to define as 2W(D t ). 
Thus the order of the Boussinesq equation (0) is 4, the order of equation JSJl is 3. 
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2.3 Recursion, Hamiltonian and symplectic operators 



Integrability of an evolutionary equation we associate with the existence of an infinite hierarchy of its symmetries. 
If a recursion operator is known, then symmetries can be found recursively and therefore the existence of a 
recursion operator is a sufficient condition for integrability. 

For evolutionary system (|14|> a recursion operator can be defined as a "pseudo-differential operator" R with 
coefficients in TZ which satisfies the following operator equation 

D a (R) = a*oR-Roa r , (31) 

where a* is a Frechet derivative of vector a = (P, Q)^ r . If action of R is well defined on a symmetry bi, 
i.e. b2 = -R(bi) £ IC, then i>2 is a new symmetry of the evolutionary system 1)1 4[l. Starting from a "seed" 
symmetry bi, one can build up an infinite hierarchy of symmetries b n +i = R n (bi), provided that each action 
of R produces an element of the ring TZ, i.e. R n (hi) £ IC. 

For example, if we represent the Boussinesq equation J2J) in the form of evolutionary system (JSJ, then 

D x 
Dl + 2uD x + 2u x 

and it is easy to verify that a pseudo-differential operator 

/ 3v + 2v 1 D- 1 ADl + 2u + u x D- 1 \ 

\ AD* + WuDl + 15uiD x + 9u 2 + Au 2 + (2u 3 + Auux)D~ 1 Sv + v^- 1 J {6 ' 

satisfies equation 1)31)1 and therefore is a recursion operator. 

In the previous section we have shown that the second component of a symmetry is completely determined by 
its first component. We can restrict the action of the recursion operator on the first component. If S is the first 
component of a symmetry of the Boussinesq equation JHJ, then 

S = R(S) = (3w + 2v 1 D- 1 + (AD 2 X + 2u + uiD~ 1 )D t D~ 1 ) (S) 

is the first component of the next symmetry in the hierarchy. We call 5ft a restricted recursion operator. The 
restricted recursion operator 5ft of the Boussinesq equation is a homogeneous pseudo-differential operator of 
weight WQR) = 3 (the weight W^P" 1 ) = -1), therefore ord(S) = ord(S) + 3. 

A space-shift, generated by bi = (ui,vi)^ T , is a symmetry of the Boussinesq equation (JSJ. Taking u% as a 
seed, we can construct an infinite hierarchy S^k+i = $t k (Sx) of symmetries of orders 3fc + 1, k = 0, 1, 2, . . .. For 
example 

S4 = 4«3 + Av\u + Avu\ = AD x {v2 + vu) . 

We see that S4 is a total derivative and therefore SV = 5ft(54) £ 1Z is the next symmetry in the hierarchy, etc. 
The Boussinesq equation itself is not a member of this hierarchy. If we take a seed symmetry, corresponding to 
the time-shift Ci = (vx, W3 + 2uui)^ T we receive another infinite hierarchy of symmetries = 5ft fc (vi), k — 

0, 1, 2, . . .. The Boussinesq equation does not have symmetries of order 3k, k E N. One can show that ^fe+i and 
S3k+2 are elements of the ring 1Z for any k £ N and therefore 5ft generates two infinite hierarchies of symmetries 
of the Boussinesq equation. Moreover, all symmetries from the both hierarchies commute with each other. 

For system a recursion operator is completely determined by its two entries i?n and R\2- Indeed, it follows 
from <|32) that 

R21 = D- 1 o (D t {Rxx) + Rx2 o G*u ) , R22 = D- 1 o (D t {R 12 ) + R n o D x + R 12 o G w ) . (33) 
The restricted recursion operator Jft for system (JHJ can be represented as 

5ft = R 11 +R 12 D t D~ 1 . (34) 

If the system of equations is bi-Hamiltonian, then there is a useful splitting of the recursion operator in a 
composition R = TL o J of a Hamiltonian H and symplectic J operators. A comprehensive definition and 
theory of these operators one can find in the fundamental monograph of Irene Dorfman |14|. Here we remind 
basic facts and introduce notations. 

A "pseudo-differential operator" H is called a Hamiltonian operator if it is skew-symmetric TV = —Ti and 
defines a Poisson bracket for any two elements /, g £ TZ: 



{f,g}=(sf tr -nsg)) 
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which satisfy the Jacobi identity. A Hamiltonian operator H of an evolutionary system (|14f> should also be 
invariant (its Lie derivative should vanish): 

D a H = a»oH|Hoa!. (35) 

An evolutionary system 1)14(1 is called Hamiltonian, if a = H.(5H), where H is a density of a conservation law. 
The density H is called the Hamiltonian of the system. If p is a density of a conservation law of a Hamiltonian 
system l|14jl. then b = TC(Sp) is a generator of a symmetry of the system. 

A symplectic operator J is, in a sense, inverse to a Hamiltonian operator. It is also skew-symmetric J> = — J 
and defines a closed 2-form. It should also be invariant (its Lie derivative should vanish): 

D a J = -at o J - J o a* . (36) 

It follows from 1(35(1 . I|36(l that the composition Ho J satisfy equation 1)31(1 . Indeed, 

D a {H oJ) = D a (H) oJ + Ho D a {J) =a,oHoJ + HoatoJ-HoaloJ-WoJoa t . 

In the case of evolutionary system (JSJ) a Hamiltonian operator 7i is 2 x 2 matrix whose entries Hij are differential 
(or pseudo-differential) operators. It follows from ((35(1 that 

H 22 = 0™ o H u o D" 1 + GU o W 2 i o Z}- 1 - £> a (W 21 ) o .D- 1 (37) 

and W12 = — H|i due to the skew-symmetry. Thus two entries Mn and H21 completely determine the Hamil- 
tonian operator. 

Similarly for a symplectic operator we find from ((36(1 that 

J11 = -GL o J 22 o C- 1 - J 12 o o i?- 1 - D a (J 12 ) o i?- 1 (38) 

and J 2 \ = —Jl 2 - Two components J 22 and J7i2 completely determine the symplectic operator of system 
A restricted recursion operator for equation 10 can be represented in the form 

3? = Tin o Jn + H12 ° J21 + (Hn ° Jia + W12 ° J22) o D x 1 D t . (39) 
In the case of the Boussinesq equation (JHJ) there are Hamiltonian and symplectic operators of the form 



3D X + uD x + D x u 3vD x + 2v ± \ _ ( D~ l 

3vD x + V! ADl + 5(uD 3 x + D 3 x u) - 3{u 2 D x + D x u 2 ) + AuD x ii J ' V f^ 1 



It is easy to check that J 1 is also a Hamiltonian operator for the Boussinesq equation and compatible with 
H. 

2.4 Symbolic Representation 

In this section we remind basic definitions and notations of perturbative symmetry approach in symbolic rep- 
resentation (for more details see e.g. [9"lll0|'). 

We start with the symbolic representation 1Z of the ring 1Z. Linear monomials u n , v m are represented by 

u n -> uff, v m -» «CT (40) 
Quadratic monomials u n Um, u n ^m, ^n^m have the following symbols 

yCOT+ffm -> en ^ m -y(crc 2 m + crc2") («) 

A symbolic representation of a monomial Uq ?/™ 1 ■ • ■ itp p v™° v™ 1 ■ ■ ■ v™ q , no + n± + • — \-n p = n, mo + mi + • • • + 
m g = m is defined as: 



— > It U (^1^2 ' ' ' £n £n + l ' ' ' £n +ni ' ' ' (C1C2 ' ' ' Cm Cm Q + l ' ' ' Croo+TOl ' ' ' Cm)c ! (42) 
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where triangular brackets (}j and ()q denote the averaging over the group of permutations of n elements 
£i> • • • >£n, denoted by £„, i.e. 

(c(£i,...,£ n ,Cl,---,Cm))f = — } ^2 c(cr(£i), . . . ,cr(£„),Cl, ■ ■ • , Cm) 
and m elements Cl, •••> Cm respectively. 

To define the symbolic representation 1Z of the ring 1Z we need to define summation, multiplication and derivation 
in 1Z. To the sum of two elements of the ring corresponds the sum of their symbols. To the product of 
two elements f,g G K with symbols / -» u"u m a(£i, . . . , £„, Ci, • • • , Cm) and g -> u p v q b{^!, . . . , f p , Ci, • • • , (q) 
corresponds: 



fg - U "+^ m +«((a(Ci, . . . , £„, Ci, ■ ■ ■ , CmMn+i, • • • , U+p, Cm+1, • • • , Cm+,)Mo (43) 

where the symmetrisation operations are taken with respect to permutations of all arguments £ and C- It is 
easy to see that representations of quadratic (|41l) and general (|42|l monomials immediately follows from (|4(J|) 
and gHJl- 

If / £ 7?. has a symbol / — > u n v m a(£i, . . . , £ n , Ci, • • • j Cm), then the symbolic representation for its N-th derivative 
D N (f) is: 

D%(f) -> uVfo + 6 + • • ■ + £„ + Ci + C 2 + • • ■ Cm) JV a(a, ...,£«,&,..., Cm) 
Thus we obtained the symbolic representation 1Z of the differential ring 1Z. 

To construct the symbolic representation JZ\rf\ of the associative ring of formal series 1Zc[D x ] we define a symbol 
•q which corresponds to the derivation operator D x and satisfies the following rules of action and composition: 

r 1 (u n v m a(Z 1 ,...,Z n ,(i,...,( m )) = 

= U n V m (^ +6 + •• • + £„ + & +C2 + • --CmMfc. ■ • ■ ■ ■ ■ ,Cm), 

7/OU n t; m a(^l,...,Cr i ,Cl ) ---,Cm) = 

= U n U ro (Cl + 6 + • • • + Cn + Cl + C2 + ■ • • Cm + T?)a(&, ■ ■ ■ , £n, Cl, • • • , Cm) 

The last expression represents the Leibnitz rule D x o f = D x (f) + / o Z) x in the symbolic representation. 
Moreover, for any iV 6 Z 

D N of ^ ou » t ,™ a ^ lj . . . ,^ n , d, . . . , Cm) = «"«" l (Cl+6 + - ■ -+£n + Cl + C2 + - ■ ■ Cm+^a^l, •-.,£„, Cl, Cm)- 



For any two terms /-D£, <7-D| of formal series (p, q € Z) with symbols / — ► w™v m a(£i, . . . , £„, Ci, ■ ■ ■ , Cm), g — * 
u s v r b(^i, . . . , £ s , Ci) • • • , Cr) the composition rule in the symbolic representation reads 

fDlogDl^ 

U n+S V m+r ((a{t U . . . , £„, Cl, ■ ■ • , Cm)(£n+1 + C-»+2 + • • • + Cn+s + Cm+1 + Cm+2 + ' • • + Cm+r + vT X 
X&(£n+l,£n+2, ■ ■ • ,£n+s, Cm+1, Cm+2, ■ • ■ , Cm+r)?/ 9 )^, (44) 

where the symmetrisation is taken with respect to permutations of arguments £ and arguments C, but not the 
argument rj. 

Elements of lZ[n] are formal series 

A = a 00 {rj) + uaw (£1 , rj) +vaoi (Ci , rj) +u 2 a 2 o (£1 , £2, rj) +uva n (€i , Ci , rj) + u 2 a 02 (Cl , C2 , f?) +u 3 a 30 (£1 , f 2 , £3, »?) H , 

(45) 

where functions a nm (£i, ...,£„, Ci, ■ ■ ■ , Cm, rj) are symmetric functions with respect to permutations of arguments 
£i and arguments Ci- Functions a nm (£i, . . . , £ n , Cl, • • • , Cm, rj) are not necessarily polynomials, but they possess 
an important property which we call locality: 

Definition 2. We shall call a function a nm (£i, Ci, ■ • ■ , Cm, rj) local if all the coefficients of its expansion 

OnmiZl, ■ ■ ■ ,£n,Cl, ■ ■ - ,Cm,7?) = 5Z a «™(^' ' ■ ' >£™>Cl, ■ ■ ■ ,Cm)?7 fe , ?7 ~> OO 

are polynomials in variables £1, • • • , £«, Cl, • • • > Cm- 
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The property of locality reflects the fact that coefficients of a formal series A £ TZ[D X ] are differential polynomials, 
i.e. elements of the ring TZc- 

Formal series (|45ll form a ring: the summation is evident, while the multiplication (composition) rule is given 
by (compare with l|44|0 : 

u n v m a nm {£, 1 , . . . ,£„,Ci, ■ ■ -Xm,v) ou p v r b pr (£, 1 , . . . ,£ p ,(l, . . .,Cr,f?) = 

= U n+p V m+r {{a nm (£i, ...,tinXl,---,(m,V + £n+l H 1" fn+p + Cm+1 H h Cm+r) X 

xfrpr(£n+l, . . . , £n+p, Cm+1, • ■ ■ j Cm+r, V)) i) C ■ 

The ring 7^ [77] is graded 

where 72-" [77] with n = 0, 1,2, etc. are constant (i.e. u, v independent), linear in u or v, quadratic, cubic, etc. In 
other words, TZ n [rj\ is eigenspace of operator D , cf. ltTT|) corresponding to eigenvalue n. We say that a formal 
series A = o{TZ n [rj\) if A £ fe>n ft fc fo]. 



2.5 Formal Recursion Operator 

Choosing symbolic representation we make calculus of derivations much simpler paying the price of a complicated 
multiplication rule for symbols. For example, the symbolic representation of the Frechet derivative /* = 
(f*,u,f*,v) of an element / £ TZ with a symbol u n v m a(£i, £ n , £i, . . . , Cm) is: 

(/*,«, /*,«) -> (ntt n_1 t) m a(Ci, . . . , £ n -i,V, Ci, • ■ ■ , Cm), mifV^a^i, . . . , £„, Ci, ■ ■ • , Cm-i, ??)) ■ (46) 
The variational derivative l|25|) of / in symbolic representation takes form 

(n— 1 m— 1 \ 

nw n - 1 W m a(^i,...,C„-l,-E^'Cl:---:Cm),TOU™7j m - 1 a(Cl,...,Cn,Cl,---,Cm-l, - E ^) I • 
fc=l fc=l / 

In symbolic representation the problem of computation of symmetries and conservation laws can be reduced to 
linear algebra [TU1I§|. 

We start with system ^ where G £ TZ and it can be represented as 

G = Gi + G 2 + G 3 + ■ ■ ■ , G n £TZ n , 

and the corresponding evolutionary derivation is generated by the vector a = (t>i,G)^ r . 
Its symbolic representation is of the form 

u t = v(i 

v t = + vuj 2 {(i) + u 2 a2o(Ci,6) + uwan(Ci,Ci) + w 2 ao2(Ci, C2) + « 3 03o(6>6>6)+ (47) 
+u 2 fa2i(Ci,6,Ci) + ™ 2 ai 2 (Ci,Ci,C2) + w 3 a 03 (Ci,C2,C3) H = G 

Symmetries of system (jHJl , which are generated by elements of TZ, can be represented by evolutionary equations 

u T =un 1 (£ 1 ) + vn 2 (Ci ) + u 2 A 20 (£1 , &) + uvA n (fr, Ci) + v 2 A 02 (Ci, C2) + u 3 A30 (Ci ,6,6)+ (48) 
+u 2 ^2i (6 , 6 , Ci ) + "« 2 ^i2 (Ci , Ci , C2) + v 3 A 03 (Ci , C2 , Ca) + • • ■ 

where f2i(£i), S^2 (Ci) , ^ij are polynomials. 

It follows from l|4tj[l that the Frechet derivative a* in the symbolic representation has the form 

a*=(/ J ) (49) 



G* r u G*^ v 
where 

G*,„ = u7i(jj) + 2ua 20 (£i,?7) + «oii(J7,Ci) + 3u 2 a 30 (6, £2, v) + 2wra 2 i(£i, r), Ci) + v 2 a 12 {ri, Ci, C2) + 
G»,„ = w 2 (?7) + wan (^1,77) + 2i;a 2 (&>*?) + « 2 a2i(Ci, 6, 17) + 2waia(£i, Ci, + 3v 2 a 03 (Ci, C2, v) + 
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Let us define a 2 x 2 matrix R whose entries are formal series 

Rll = <Poo(v) + «</>io(£i,»7) + -v0oi(Ci,»7) + m 2 02o(Ci,6,?/) + W^n(£i, Ci, ??) + " 2 0O2(Ci,C2,'7) H (50) 

Hi 2 = Voo(»7) + 2 U ^io(^i)??) + WtoiCCi* 7 /) +u 2 ip2o(Ci,£,2,v) +uvipxx((i,Ci,r]) + v 2 ip 2((i, ( 2 , r{) H 

and (compare with 

i?21 =^ 1 o(i?n, f +A 12 oG,, u ) ) i?22 = O (^12,4 + i?12 O G*,„ + £ U O 7?) (52) 

Definition 3. We shall call matrix R a formal recursion operator of system \4-l\j if it satisfies equation 

Rt = [&*,R] (53) 
and all the coefficients 4>ij , ipij > h 3 = 0) 1> 2) • • • of the formal series R\\ and R\% are local. 

Any recursion operator in symbolic representation satisfies the above Definition. The category of formal series 
is less restrictive than the category of pseudo-differential operators, since the action of formal series on the ring 
1Z is not defined and there is no issue of convergence. Actually we always consider a few first terms of i?n and 
R\2- Therefore it is convenient to introduce a notion of an approximate formal recursion operator. 

Definition 4. We call matrix R an approximate formal recursion operator of degree n if it satisfies 

R t -[a„R}^o(TZ n [ V }) 
and all the coefficients up to degree n of formal series R\\ and R\2 are local. 

Theorem 1. Suppose that system |^7| ) possess infinitely many approximate symmetries \4$il of degree m 

u Ti = ufii,i(£i) + wfi 2 ,i(Ci) + o^ 1 ) • 

• If < deg(f7i j i) < deg(f2i.2) < deg^i^) < . . . , then system |^7| ) possess an approximate formal 
recursion operators of degree m — 1 with (poo{f]) = Vi ^oo(^) = 0. 

• If < deg(f22,i) < deg(f22,2) < deg(f22,3) < ... , then system J^7] ) possess an approximate formal 
recursion operators of degree m — 1 with 0oo(?7) = 0, ipoo(n) = r). 

The proof of the Theorem is straightforward and very similar to the proof of the existence of a formal recursion 
operator given in UJ. 

In symbolic representation equation (|53|l can be solved and all the coefficients 4>ij,ipij, i, j = 0, 1, 2, 3, . . . can be 
expressed in terms of the right-hand side of system 147JI . i.e. u>i, a{j. A substitution of 149(1 . H5()[1 . H51|1 . I|52() in l|53(l 

leads to a system of linear algebraic equations on the coefficient functions of the formal recursion operator. There 
are no restrictions on the choice of (/>oo(??), ipooiv) an d they can be chosen arbitrary (for example as in the above 
Theorem: O o(»?) = V, 4>oo(.v) = or (f>o (v) = 0, tpoo(v) = v)- To determine coefficients <f>ij(€,'n),ipij(€,'n) 
when i + j = 1, we have to solve a system of four linear equations. 

In what follows we will study homogeneous equations Q of even order p = 2n, which can be written in the 
form ©. Thus in (|47Jl we assume 

where a, are constant parameters. Without loss of generality we can choose 

If p, 7^ 0, ±/3 then the solution of equations for </>io(£> if), ipia(£, rj), </>oi {£,, rj), i>oi(£, rj) can be written in the form 
08 - M) 2 /i(M, t v) + W + ^) 2 /i(-M, t, V) + (^ 2 ~ /3 2 )(/2(^ ; e, »7) + / 2 (-/i, V)) , 

([3 - (m, e, 77) + (/3 + 6 »?) + (/*- /3)/ 2 (m, £, - (M + /?)/ 2 (-M, £, »7)) , 

(/? - m)/!^, e, v) + (P+ n)fi(-n, & »?)-(/*+ /?)/2(m, e, »7) + - /w-M, e, ^)) , 



<l>oi(£,v) 



2<Tn-l 



^01 (^,7?) = 



1 



4^2^-l^r, 
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where 
and 



52 (a*) £, v)d3(—^, & v) 



<S + v) 



2(<M£ + V)~ Mv)) + + t*)((v n + C)1>oo(S + v)~ V n Tpoo(ri)) 



'00 



(£ + v) - Mv)) + M)*7 n (V>oo(£ + »?) - V-ooW) + + »)ei>o (Z + V) 



gxi^v) = (P-m+vr-(P+m n +V n ), 

g 2 (^tv) = W-vM + fl) n -(f3 + riC-(P~riv n , 
9M,V) = W-vM + v) n -(P-riC-(P + riri n , 

gt(n,z,v) = iP-mi+v) n -e-v n ) 

n (/x, f , n) = 4a 20 (£, 77) + 03 + /^""Wfc »?) + (/? + /*)f -1 on fa, + {fi + A*) 2 ? 1- V^M*, V), 
r 2 (/x,e,r?) = 4a 20 (e,»7) + (/?- mX" 1 an + ifi + /^"Vifa.f) + (/3 2 - ^""V" 1 ^^) 

In a similar way, solving a linear system of 8 equations, one can find coefficients 4>20, 4>iij 4*02, V^O) ipu, V'02, etc. 
Despite of a straightforward way of computation, the explicit expressions for the coefficients ^20 , ■ ■ ■ 1P02 are 
quite cumbersome and, therefore we do not write them down in the paper. 

Cases a — (i.e. /i = ±/3) and 4a + 1 = (i.e. /i = 0) require a consideration of higher degree terms and have 
to be treated separately. 



3 Integrable equations of 4th, 6th and 10th order 

In this section we consider homogeneous partial differential equations of the form 

u t =vi, v t = cm 2n -i + f3v n + G(u, V, . . .U 2n -2, v n -i), (54) 

with n = 2, 3, 5, assuming that the polynomial G(it, v, . . . U% n —2, Vn—i) has non-zero quadratic terms and param- 
eters a and (3 do not vanish simultaneously. We also assume that the weight W(u) of the variable u is positive 
(it is automatically integer, since G has a quadratic part). It follows from 1)54(1 that the weight of the variable 
v is W(v) = W(u) +n — 1. In the case dG/dv — 0, system Q54JI corresponds to a non-evolutionary equation of 
form Q with of order p = 2n (and q = n). 

We select equations <|54ll that satisfy necessary conditions for the existence of higher symmetries, i.e. the 
conditions that a first few coefficients of the corresponding formal recursion operator are local functions and 
then prove their integrability (the existence of an infinite hierarchy of higher symmetries). For all equations 
found we show that they satisfy sufficient conditions for integrability, such as they possess bi-Hamiltonian 
structures, recursion operators and the Lax representations or linearising substitutions. 

3.1 The 4th order equations 

It is easy to see that a homogeneous system (|54f) with n — 2 is linear if W(u) > 3. In the case W(u) — 3 the 
only possibility is G = 7U 2 , 7 / which leads to a non-integrable equation for any choice of a, j3 and 7. Thus 
the the weight W(u) can be equal to 2 or 1. 

3.1.1 The case of W(u) = 2 

The most general nonlinear homogeneous system of equations 1)54|1 with W(u) 
is of the form: 

f u t =vi 

Vt — au 3 + I3v 2 + c\UU\ + c 2 uv, 

where ci, c 2 are arbitrary constants and at least one of them is not zero. 



= 2 (correspondingly W(v) = 3) 

(55) 
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Without loss of generality we need to consider the following three types of system 

u t = vi 



2 — ft 2 

v t = M 4 u 3 + /3t?2 + ciuui + c 2 uv, u {0, ±/3}, ^,/3eC, 



(56) 

J U-3 T ^1^2 T t-iu.u.1 T C2U-U) ft f \ u j - s -h' il ftF c 

Ut = Wl + , (57) 

Vt = V2 + C\UUi + C2UV 

(58) 

V t = -4U3 + V 2 + ClUUl + c 2 uv 

The first system represents the generic case a £ {0, — ^-}, while the other two represent the degenerate cases. 

Proposition 1. System f,56)) possess two formal recursion operators with 0oo(??) = "7 ipooiv) = and (j>Qo(v) = 
0, tpoo(rj) — t) if and only if = C2 = 0. -By re-scalings it can be put in the form 

Ut=V \ 9 (59) 

Systems \5l\) and $58p are not integrable, they do not possess a formal recursion operator with <^>oo = 0, V'oo = V 
unless Ci = C2 = 0. 

System ()59|l represents the Boussinesq equation (J2J, which is known to be integrable. Its Lax representation 
can be found in Recursion, Hamiltonian and symplectic operators for the Boussinesq equation were given 
in Section 12.31 In the Proposition 1 and below by re-scalings we mean an invertible change of variable of the 
form: 

u — > ol\ u, v— >a 2 v, x^-a^x, t— > cut, ati G C . 
3.1.2 The case of W(u) = 1 

The most general homogeneous system of equations (15411 corresponding to the case W(u) = 1 and n = 2 is: 

ut = vi , . 

v t — auz + (3v 2 + ciuu 2 + c 2 u\ + C3U1V + C4UV1 + c^v 2 + cqu 2 u\ + c-;u 2 v + c$u 4 

where Ci, i = 1 ... 8 are arbitrary constants and we assume that at least one of the coefficients ci, . . . , C5 is not 
zero. Without loss of generality we consider the following three types of the system (|6Ull : 



u t = Vl 



Vt - 4 



M 4^ U3 + /3«2 + C1MM2 + C 2 Wi + C3U1V + C4UV1 + C 5 V 2 + CqU 2 Ui + C7U 2 V + C 8 U 4 , 

where \i ^ {0, ±/3}, fx, (5 G C and 

M t = Ul 

ft = 1' 2 + C1UU2 + C 2 U 2 + C3U1V + C4UV1 + C5V 2 + CqU 2 Ui + C7U 2 V + CsU 4 
u t = Vl 

V t — ~jU 3 + V 2 + C1UU2 + C 2 u\ + C 3 UiV + C4UV1 + C 5 V 2 + C 6 U 2 Ui + C 7 U 2 V + C&U 4 



(61) 

(62) 
(63) 



Proposition 2. System possess two formal recursion operators with (poo(T]) = rj, - 0oo( ? 7) = and 0oo( 7 ?) = 
0, iftooiv) =r l if an d on ly */ ( U P t° re-scalings) it is one of the list 

Ut = Vl ,2 (64) 



v t = u 3 + u\ 

U t = Vl 

vt = U3 + 2uiv + 2u 2 ui 
u t = Vl 

vt = us + 2uiv + 4uvi — 6u 2 u 



a... ( 65 ) 



a... ( 66 ) 



u t = Vl 

vt = u$ + 4uu 2 + 3u 2 — v 2 + 6u 2 ui + u 



2 ,;2 r , ,,-> ' (67) 



U t = Vi ^ , . 

v t = aus + V2 + 4auu2 + 3au 2 + mv + 2uvi — v 2 + 6au 2 ui + u 2 v + au 4 , a ^ —\ 
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System possess a formal recursion operator with (ftooiv) = 0> Tpooiv) = V if an d on ^V if ( U P t° re-scalings) it 
is one of the list 

' Ut = V \, (69) 

Vt = V2 + LUV\ 



u t = Vi 

V t — V2 — u\ + 2u\V — V 



(70) 



u t = V! 

vt = V2 — 2uu2 — 2u\ + 2u\v + 6uvi — \2u 2 u\ 

System I6',9j) possess two formal recursion operators with (f>oo(r]) = rj, ipooii]) = and (pooiv) = 0, i>oo(r]) = r] if 
and only if (up to re-scalings) it is 

Ut = Vl (72) 
Vt = + V2 — UU2 — \u\ + u±v + 2uv\ — v 2 — |m 2 ui + u 2 v — \u A . 

Only equations 164(1 and l|69|l from this list represent second order equations of the form in the standard set 
of dynamical variables, since their right hand side do not depend on the variable v explicitly. Equation (|f>4|) 
is a potential version of the Boussinesq equation, cf. (|28|l . Indeed, if we introduce variables U\ = u, V\ = v 
then in terms of U,V equation (|59|l takes form I|tj4fl . Its symmetries, conservation laws, recursion and multi- 
Hamiltonian structures can be easily recovered from the theory of the Boussinesq equation. For example, the 
recursion operator of the potential Boussinesq equation (|64ll is D^ 1 o Ro D x where R is the recursion operator 
(|32|l in which Uk,Vk is replaced by Uk+i,Vk+i and therefore the corresponding restricted recursion operator for 
equation 164J1 is of the form 

= -DjT 1 o (3v-lD x + 2v 2 + {AD 2 X + 2u x + u 2 D~ l ) o D t ) . 

Equations (|65|l and l|66() are known to be integrable. Corresponding Lax representations and references can be 
found in |17j . 

Equations (|67[) and l(68|l can be mapped into linear equations 

w tt = Wi , w tt = awi + W2t 

respectively by the Cole-Hopf transformation u — (logiy) x . In these cases the restricted recursion operator is 
the same as for the Burgers equation 

3? = D x + u + u^D- 1 , (73) 
which generates higher symmetries from the seeds u\ and u t . 

Equation l|69(l can be reduced to the Burgers equation with time independent forcing 

u t — U2 + 2uu\ + wx , wt — 
by invertible transformation v = w + u\ + u 2 . The latter can be linearised by the Cole-Hopf transformation. 
Equation l|7U|) can be reduced to the system 

Ut — V.2 + U!l , Wt = — W 2 

by a simple invertible the change of the variable v = w + u%. System l(7()|l provides an example of an equation 
that possesses neither higher symmetries nor a recursion operator. However, a formal recursion operator does 
exist and therefore it is in the list of the Proposition. Its integration can be reduced to the integration of a 
linear nonhomogeneous heat equation with a source term of a special form. 

System (|71|) possesses an infinite hierarchy of symmetries of all orders generated by a recursion operator 

$ = -2u + 2u 1 D- 1 +D t D- 1 

starting from the seed u\. 

By a simple shift of the variable v — w + u\ + 2u 2 , system (|71f) can be transformed in the form 

Ut = U2 + 4uui +idi, Wt = 2D x (uw) . (74) 

This system can be linearised by a reciprocal and then point transformations (we would like to thank A. Hone 
who helps us to find an explicit form of the linearising transformation). 

Equation lj721l is a particular case of corresponding to the exceptional case a = — \ . 
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3.2 The 6th order equations 



Homogeneous 6th order (n = 3) equations l|54|) with non-zero quadratic terms correspond to W(u) < 5. We 
restrict ourself with the case W(u) > 0. The weights 3, 4, 5 do not lead to integrable equations: 

Proposition 3. For weights W{u) = 3,4,5 there are no equations possessing a formal recursion operator. 
3.2.1 The case of W(u) = 2 

The most general homogeneous system 1)540 corresponding to W{u) — 2 can be written as 

u t = vi (?5 , 
v t — au 5 + f3v 3 + D x [ciuu 2 + c 2 ul + c 5 u 3 ] + c 3 uvi + C4VU1, 

where ct,f3, Cj, i = 1, . . . , 4 are arbitrary constants and we assume that at least one of ci, . . . , C4 is not zero. 
Without loss of generality we consider the following three types of the above system: 

Ut=V \ 2 (76) 
v t = M ^ u 5 + /3v 3 + D x [ciuu2 + c 2 u\ + c 5 u 3 ] + c 3 uvi + c^vui, /i^{0,±/3}, /i,/3eC 



Vt = V 3 + D x [c\UU 2 + C 2 u\ + C 5 U 3 ] + C3UV1 + C4VU1 



Ut = Vl (78) 
v t = —JU5 + V3 + D x [c\UU2 + c 2 uf + c 5 u 3 ] + c 3 uv\ + C4VU1 

Proposition 4. If system |76p possess two formal recursion operators with </>oo(?7) = ij, Voo(?7) — and (j>oo(rj) 
0. tpoo(j]) = r\ then, up to re-scalings, it is one of the list 



u t = Vl 



Vt = 2u 5 + v 3 + D x [2uu 2 



u\ + ±u 3 ] 



(79) 



Ut = Vl (80) 
Vt — ^u 5 + v 3 + D x [uu 2 + uv + iy 3 ] 

u t = v x 

v t = §1*5 + v 3 + D x [2uu 2 + \u\ + 2uv + §w 3 ] v > 

If system {77}) possess a formal recursion operator with (fiooiv) = 0, ^Pooiv) — V then, up to re-scalings, it is one 
of the list 



u t = Vl 

vt = v 3 + uvi + u\v 
u t = Vi 

vt = v 3 + 2uu 3 + Au\u 2 — Au\v — &uv\ — 2Au 2 u\ 



(82) 
(83) 



System \78ty does not possess a formal recursion operator for any non-trivial choice of Ci 
Equation 1)79)) can be rewritten in the form Q): 



u tt = 2u 6 + u 3jt + D 2 X ( 2uu 2 + u\ + ^-u 3 
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It has the following restricted recursion operator 



J? = Dt + ^Dt + ^utD-' + ^uDl + ^uDtD- 1 

13 7 4 1 

H ui D T H — u 2 H — v H — Mi D t Dl 2 . 

9 9 9 9 

Applying powers of 5ft to the seed symmetries u\ and w t , one can build up two infinite hierarchies of commuting 
symmetries generated by £4,1+1 = and Si n+3 — 5ft™(wt) respectively. There are no even order symmetries 

for equation 1)79)1 . Its bi-Hamiltonian structure is given by a Hamiltonian operator 

Hu - -D 3 x + ^(uD x + D x u) 
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Wu = \D b x + -vD x + i«i + -ul>2 + + Z U2jD;c 

and a symplectic operator 

•J12 = D x , J22 = . 

A simple change of the variable v = w — u 2 — u 2 /3 brings system (|79|) in the form 

2 2 
u t = -u 3 - -imi + toi , to t = 2w 3 + -mwi . 

Thus system ()79|l has a reduction (w = 0, i.e. v = —u 2 — u 2 /3) to the Korteweg de Vries equation. 

This system has the following Lax representation: 

2 2 5 1 1 

L = D$ + -uDl + - Ul D x + -u 2 + —u*--v, 

2 1 
A = AD% + -uD x + • 

The above mentioned reduction of the system corresponds to the case when L is the square of a second order 
operator (i.e. the square of the Lax operator for the KdV equation): 

System (|79f) and its Lax representation can be transformed to one of the cases studied in \7\. 

Recursion operators and bi-Hamiltonian structure as well as the Lax representations for the potential versions 
of equations 1)80(1 and H81|l will be given in the next section. 

Equation (|82|l is known to be integrable, its Lax representation can be found in [H]. Its bi-Hamiltonian structure 
is given by a Hamiltonian operator 

H n = 5D 3 X + 2uD x + u x 
Wia = \{Pi + uDl + 2 Ul D 2 x + u 2 D x ) + AvD x + 3« x 

and a symplectic operator 

J21 = D X + uD x \ J 22 = 2D' 1 . 



Its restricted recursion operator is of order 6. The hierarchies of symmetries can be generated from the seeds 
S\ = ui, S3 = Ut and S5 = 115 — 10«3 + D x (5uu 2 — lOuv + |u 3 )- 

Equation l|83() has the following restricted recursion operator 

Sft = -D x D t + uD 2 x + uiD x + u 2 - 10u 2 - 2v - 2uD x 1 D t - viD x Y - AuiD^u - uiD x 2 D t , 

which is of weight 4. There are two infinite hierarchies of symmetries S^k+i = 3? fe (iti) and 6*4^+3 = and 
no symmetries of even order. 

Its bi-Hamiltonian structure is given by a Hamiltonian operator 

Hii = -~Dl + 2uD x + ui 

Hi2 = ~\ D l + $uDl + 10uiD 2 x + 7u 2 D x + 2it 3 - 6u 2 D x - 12uui + 2vD x + vi 



and a symplectic operator 



J11 — D x , J 22 — . 



(84) 



If we introduce a new variable v = w + u 2 — 3w 2 , equation (|83|l can be rewritten in the form 2 

ut = U3 — Quui + wi 
w t — —2uwi — 4uiw 

It admits an obvious reduction (w = 0) to the KdV equation zt t = vi — 1*3 — 6uui. 

2 When our paper was completed we discovered that integrable system 1841 has recently been studied by A.B.Shabat he 
also has found a Lax representation for this system. 
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3.2.2 The case of W{u) = 1 



Homogeneous systems of equations (|54() with W(u) — 1 can be written in the form: 
u t = Vl 

v t — au 5 + (3v 3 + C\U2 + c 2 uiu 3 + C3UU4 + C4U2V + c 5 uivi + c 6 uv 2 + c 7 v 2 + c 8 uf + , (85) 

+C9UU1U2 + CiqU 2 U 3 + CnU 2 Vi + C12UU1V + + Ci 3 U 2 U 2 + Ci4U 3 U2 + CisU 3 V + C\qU A U\ + Ci 7 U 6 

where a, /3 € C, all Cj, i = 1, . . . , 17 are arbitrary constants and at least one of c\, . . . c 7 is not zero. 
We need to consider the following cases of the system 

u t = Vi 

v t = M ~± u 5 + (3v 3 + c\u\ + C2U1U3 + C3UU4 + C4M2U + c 5 uiVi + c e uv 2 + c 7 v 2 + c s uf+ (86) 

+C9MM1M2 + CiqU 2 U 3 + Ci\U 2 V\ + C\2UU\V + +Ci%U 2 u\ + Ci4U 3 U2 + Ci^U 3 V + C\qU U\ + C17U 6 

where \i {0, ±/3}, fi, (3 G. C, and 
u t = Vl 

vt — v 3 + c\u\ + c 2 uxu 3 + c 3 uu 4 + c 4 u 2 v + C5U1W1 + c 6 uv 2 + c 7 v 2 + c%u\+ (87) 

+CgUUiU 2 + CiqU 2 U 3 + C\\U 2 V\ + C\2UU\V + + Ci 3 U 2 U 2 + Ci4U 3 U2 + Ci$U 3 V + C\qU U\ + Ci 7 U 6 
U t = Vl 

v t = —\u 5 + V3 + c\u\ + c 2 uiu 3 + C3UU4 + C4U2V + + c 6 uv 2 + c 7 v 2 + c%u\+ (88) 

+C9UU1U2 + CiqU 2 U 3 + C\\U 2 V\ + C\2UU\V + + Ci 3 U 2 U 2 + Ci4U 3 U2 + Ci^U 3 V + C\qU A Ui + C17M 6 

Proposition 5. // system J^Hh}) possess two formal recursion operators with (f>oo(n) = Vi ipooiv) — 0j an d 
(f'ooiv) — 0: V'oo(^) = U P t° re-scalings, it is one of the equations in the following list 

(89) 
(90) 
(91) 



Ut 


= Vl 










v t 


= 2u 5 - 


f v 3 + 


u\ + 2u ± u 3 + ^u\ 








= Vl 










Vt 


= l U 5 


+ v 3 4 


Uiu 3 + uivi + iitf 






u t 


= Vl 










Vt 


= \u b 


+ V 3 + 


■ 2mu 3 + \u\ + 2uivi + ~u\ 






u t 


= Vl 










Vt 


= au 5 + v 3 -{ 


- 10au 2 + 15auiu 3 + 6auu4 + v%i2 + 3uiVi + 3uv2 


— v 2 + lhau\ + 


15au 2 u 3 + 




+ 60auuiU2 + Suuiv + 3u 2 vi + A5au 2 u 2 + 20au 3 U2 + u 3 v + 


15au 4 ui + au 6 , 




u t 


= Vl 










Vt 


= "5 4 


- 6UU4 


+ 15miW3 + IOm 2 — v 2 + I5u 2 u 3 + I5u 3 + 6OUU1U2 


+ 4:5u 2 ul+ 






+ 20u 3 u 2 - 


h 15u 4 ui + u 6 







(93) 

// system \8T\j possess a formal recursion operator with 4>oo(v) = 0? ^ooiv) — V> U P t° re-scalings, it is one of 
the list 

' Ut = V \ (94) 

V t = V 3 + UiVi 



U t = Vl 

vt — v 3 + SuiVi + 3uf2 + 3u z vi 



2.. (95) 



u t = Vl 

Vt — V 3 — u\ + 2U2V — V 



2^o_,_.,2 (96) 



If system $88p possess two formal recursion operators with 4>oo(v) — 7 h ipooiv) — 0j an d ^ooiv) — 0, ipooiv) = 
up to re-scalings, then it is 

u t = Vl 

An.. _L _ 5„,2 _ 15„,.„,„ _ 3„„ . , .,„,_ , o„,„, , q„,„,„ _ „2 _ 15„3 _ 15, ,2 



v t = -4U5 + v 3 — |t*2 — j-uiu 3 - |uu 4 + vu 2 + 3uiVi + 3uv 2 — v — fuf — j-u 2 u 3 + (97) 

15UUIU2 + SuUiV + 3u 2 Vi — ^-U 2 U 2 — 5u 3 U2 + U 3 V — ^j-U 4 Ui — ju e 



Equations (|89I91^ and (§Ul are "potential" versions of (|79I81^ and |jH2J 

System H9()|l represents in the form of equation as 

1 „ / 1 



u u = |t*6 + u 3t + D x I mu 3 + uiu t + |tti ) (98) 
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It is a bi-Hamiltonian system with a Hamiltonian operator 

Hn = 6D x + 2u 1 D x 1 + 2D- 1 o Ul 
33 

W21 = —D 3 x + llu 1 D x + 2v 1 D- 1 +6D- 1 ov ll 
o 

and a symplectic operator is defined by 
3 

J12 = -U!Dl + 3u 2 D* + [5u 3 + ul + 5v!]Dl + [2u 4 + 2u 1 u 2 + 5v 2 ]Dl + 
5 

+[— 2it 5 - 2mim 3 — 2mo + 4u 3 + 3uiVx]D x - -u 6 - 2u 2 u 3 - uiu 4 + v 4 + it 2 vi + ^1^2, 

5 

27 11 
J22 = -3- J D^ + 3uiZ)3+£>3o3u l + [-2u3 + -u?+ui]-Dx+-D !C °[-2u3+2 U i + v i]- 

A restricted recursion operator for equation l|98|l can be expressed in terms of H and J (see l|39t ). This system 
does not have symmetries of even order and of order 5 mod 10. The hierarchy of symmetries can be generated 
by the restricted recursion operator acting on the seed symmetries S_i = 1, Si = ui, S3 = Ut and 



14 35 35 2 70 35 35 35 

-^UlU b + —U 2 U A + — M3 + — u 3 u t + —u 2 u u + —uiu 2t + — 1 



35 9 35 35 2 

A Lax representation of system (|9"0j) is 

L = D 5 x + -u 1 D 3 x + -u 2 D 2 x +[-^u 3 + -u 2 1 --v 1 \D x , 
A = L»^+miD x . 

System (|91|) represents in the form of equation as 

1 ^ 3 2 „ 4 „ 

w« = -pu 6 + u 3 t + D x 2uiu 3 + -u t + 2uxu t + -u 1 
5 \ 2 3 

This is a bi-Hamiltonian system with a Hamiltonian operator 

Hn = '3D X + u x D x l + D x 1 m , 
n 12 = ^-D 3 X + 4D xUl + + D^ut . 



and a symplectic operator 

9 42 21 T 1 25 

J12 = -—D 7 x -— Ul D 5 x - —u 2 D% -[-u 3 + Uuj + U Vl ]D x - [-U4 + 25uiu 2 + y u 2 ]£>* + 

+ [^ M 5 + - 2«2 - — ^«3 - 24uiUi]i) x + + 4lt 2 W3 + 2uitt 4 -V4 — 8u 2 Vx - 8uiV 2 , 

A A o 

J22 = —Dl-^iutDl + D^ + i^us-Aul-ivt^ + D^us-Aul-Av!] . 

A restricted recursion operator can be expressed in terms of Ti. and J (see (|39f) ) . It is of 10-th order. The 
system does not have symmetries of even order and of order 5 mod 10. The hierarchy of symmetries can be 
generated by the restricted recursion operator acting on the seed symmetries S_i = 1, Si = ui, S3 = Ut and 

28 or 455 2 140 or 70 70 2 

S 7 = u 7 + 7u 4t + — u x u b + 35u 2 W4 + -y u 3 + —u 3 u t + 35u 2 u lt + —uiu 2t + —u t + 

140 2 2 140 2 

H — — M1U3 + 3buxu 2 + — u-yUt ■ 

The Lax representation for equation 1|) is : 

r t,5 10 ^3 r ^2 r 35 20 2 10 , „ 10 20 5 

L = D x + Y UlD * + ™z D x + iy w 3 + yui - yViJAs + y"4 + y«i"2 - ^2, 
A = D 3 X + 2u x D x + u 2 . 
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Equations (|92[1 and (|93|l can be mapped into linear equations 

w t t = aw 6 + w a>t , w tt = w 6 

respectively by the Cole-Hopf transformation u — (loguj) x . In these cases the restricted recursion operator is 
the same as for the Burgers equation, c.f. (|73[l . 

System l|95|) represents equation of the form 

u t t = D t (u 3 + 3uu 2 + 3u\ + 3u 2 ui), 

from which follows, that Ut = u 3 + '3uu 2 + 3u 2 + 3u 2 u\ + f(x) , where f{x) is an arbitrary function. This equation 
can be linearised by the Cole-Hopf transformation u = (logiu)^. 

Equation l|96(l has similar property as system (|70|l . It can also be reduced to a triangular system 

Ut = U3 + wi, w t = -w 2 

in variables u and w — v — u 2 . 

System 1)9 7|) is a particular case of 1)92(1 in the exceptional case a = — j. 



3.3 10th order equations 

In this section we present three examples of 10th order integrable non-evolutionary equations together with 
their bi-Hamiltonian structures and recursion operators. These equations are new to the best of our knowledge. 

Proposition 6. The following systems possess infinite hierarchies of higher symmetries: 
u t = Vl 

v t = ^u 9 + v 5 + D x {3uu§ + 9uiu 5 + ^u 2 -u 4 + ^j-u 2 + 2uiVi + 4w 2 + 20u 2 w 4 + 80uuiu 3 + (99) 
+60^?, + 88u 2 u 2 + ^u 3 u 2 + ^fu 2 u 2 + ^u 5 ) 



u t = vi 

Vt = ~BZ W 9 + v 5 + ^u 7 ui + |m 6 U2 + |w 5 u 3 + j|m 2 - 5u 3 vi - ^u 2 v 2 - I0uiv 3 (100) 

45 2 75,.,.,. 75 „2,. 75,.2,. 1 45,,2,. , 225,. ,,3 1 675,,2,,2 405, ,5 
— — U§U 1 Y u l u ^ u i 4~ U 3 U 1 J~ U 2 U 3 T ~2 U l v l n ~ U 3 U 1 -\ $-U 2 U 1 iq~ u i 



Ut = Vi 

vt — + 2u 2 u§ + 6U3U4 — 6u 3 v — 22u 2 v\ — 30uiv 2 — 20uv 3 + 96uuiv + 96u 2 vi (101) 
-2D x [8u 2 u 4 + 32u Ul u 3 + I3u 2 u 2 + 24uu 2 } + 120D x [4u 3 u 2 + 6u 2 u 2 ] - 3840w 4 ui 

System 1)99)1 can be rewritten in the form Q as 

9 / 65 35 

u tt = ~^ u w + u 5:< + Ac ( 3uu 6 + 9uiu 5 + —u 2 Ui + -^u 2 + 2u\u t + 4uui :t + 20u 2 u 4 + 80umu 3 + 

„„ 9 „„ 2 256 o 384 9 2 1024 * 
+60uu^ + 88u(u 2 + —u A u 2 + —u 2 u{ + — — u 5 
5 5 125 

System 1)99)1 is bi-Hamiltonian with Hamiltonian operator: 

7 1 

Hn = — D 3 r H (uD x + D x u) 

88 * 55 v ' 



11 , 9 , , 11 , 6 2 99 1 

H12 = T7^D 7 X + —DluD 2 x + —D 2 x u 2 D x + —vD x + -D 2 x u 2 D, 



Ut 



128 x 20 x 1 40 1 55 ' 5 ' ' 1 1 

and a symplectic operator 

^2i = D x , J 22 = . 
Its restricted recursion operator is of order 6: 

9 a 21 , 12 , 17 , 14 , , 9 72 51 

= 128 x + 20 * + T Ul x + ~8 U2 x + T x + 4 U3 x + Y UUl x + 40 W4 

42 9 6 7 2,1 , , 

+8uu 2 + —u\ + -v + -D x D t + -uD~ l D t + -u 1 D x 2 D t + u t D x \ 
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Hierarchies of symmetries are of order SWi+i and Se n +5 generated from the seeds u\ and u t respectively. There 
are no symmetries of even order and of order 6n + 3, n € N. 

The Lax representation for the equation (|99|l is given by: 



32 ^4 ,„ ^3 , 96 256 2^2 ,64 768 752 256 

L = D 6 + —uD i + 16u 1 D 3 + (—u 2 + —u 2 )D 2 + {—u 3 + —uu 1 )D + —u 4 + 



.2 



2816 8192 o 256 , 16 256 256 4096 , 128 , 

H UU9 H u H w + ( Us H u-iUi H mi? H u u-t H Ui )L> 

225 3375 225 V 225 5 225 225 3 1125 225 ; 

9 32 64 

A = --D 5 - 6uL» 3 - 12 Ul D 2 - (10u 2 + — u 2 )D - 4u 3 - — uu y . 
8 5 5 

After a simple invertible change of variables 

1 „ 32 
— i 
15 

equation 1(99(1 can be rewritten in the form 

Ut — —h u 5 ~ 2uu 3 — 2uiU2 — ^ru 2 ui + Wi 

Q 32 9 

w i = 1^5 + 6ltU>3 + 6M1W2 + ^U 2 Wi + W\ 

It admits a reduction w = to the Sawada-Kotera equation ^H] ■ 

Much less is known about new integrable system l|100|) . It can be rewritten in the form (JIJ 

1 /5 5 5 25 15 

u t t = -77«io + u 5 t + D x -u 7 ui + -u 6 u 2 + -u 5 u 3 + —ul - 5u 3 u t - —u 2 u lt - I0uiu 2t 
54 V D 3 2 12 I 



v = w — —114 — 2uu 2 -it 3 , 



45 2 75 75 2 75 2 45 , 225 , 675 , , 405 

' Y U5M i ~ Y UlU2Ui ~ X"" 3 "" 1 ~~ I" 2 " 3 + Y u i u * + ""3^! + —u 2 u 1 - —% 

System 1(100(1 is a bi-Hamiltonian system. We have found the following Hamiltonian operator 
H u = ~D x + uxD~ l + D- X u x 

Hu = ~^DI + 16u 1 DI + 32u 2 DI + [yu 3 -^u 2 1 ]D x + yU4-75u 1 u 2 + 5v 1 D- 1 +D- 1 v 1 
and a symplectic operator 

Ju = -\ D l + 2u i D l + ^D 4 X + [u 3 - Y u i] D l - + [~7«5 + 30mim 3 + + 6vi ^ Dx ~ 

135 

— 3«6 + 6OM2W3 + 30ui?/4 — u\u 2 + 3^2, 

J22 = -2D 3 X + 3 Ul D x + 3D xUl 

It follows from 139(1 that the restricted recursion operator for equation 1)100(1 constructed from the above Hamil- 
tonian and symplectic operators is of order 12. The hierarchies of symmetries can be generated from the seed 
symmetries 5_i = 1, Sj. = Ui, S5 — Ut and 

105 2 ,„„ 945 2 2 835 4 

67 = U7 — 42u2t — 42uiU5 — u 3 + 189uiUt H — — u 1 u 3 — u 1 

2 2 8 

Equation 1(100(1 does not have symmetries of order 6n + 3 and even order symmetries. We have not found the 
Lax representation for this equation either. 

System ((101(1 does not represent a single non-evolutionary equation in the standard set of dynamical variables, 
since the right hand side of the second equation depends on the variable v explicitly. It can be written in the 
form Q after the extension of the set by a potential z, such that z% = u. 

System 1(101(1 is bi-Hamiltonian. We have found the following Hamiltonian 
H n =--Dl + uD x + D x u 

Hia = + 12D 5 x u - 24D xUl + 26D x u 2 - 80D 3 x u 2 - UD 2 x u 3 + l60D 2 x u Ul + 3D x u 4 

-60D x uu 2 - 45D x ul + 80D x u 3 + 3D x v - m 
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and a symplectic operators 



J21 = D. 



-1 



J22 = . 



X 



The corresponding restricted recursion operator is of order 6: 




t, 



which generates the higher symmetries from the seeds u\ and Ut- 

This equation has no symmetries of even order and of order 6n + 3. Its hierarchies of symmetries can be 
generated from the seed symmetries S\ = u\ and S5 = u t . 



Condition w = is an obvious reduction of this system. The reduced equation is the well known Kaup- 
Kupershmidt equation |20| . 

Summary 

Most of results obtained in the theory of integrable equations concern the case of evolutionary equations. 
The theory of nonevolutionary equations has new features. In order to study a nonevolutionary equation 
we represent it by a system of evolutionary equations. Such representation is not unique. A choice of the 
representation determines a set of dynamical variables and the corresponding differential ring. The concept of 
local symmetry, which is central for the theory of integrable equations, can be rigorously defined in the terms of 
this differential ring. Existence or absence of infinite hierarchies of local symmetries may depend on the choice 
of the representation, i.e. on the ring or its extension. In contrast to the evolutionary case, the existence of one 
higher symmetry does not guarantee the existence of an infinite hierarchy of symmetries |21| . 

We have developed a perturbative symmetry approach suitable for testing for integrability (testing necessary 
conditions for the existence of higher local symmetries) . It is based on the elements of the standard symmetry 
approach and the symbolic representation of the differential ring. Our approach proved to be very effective in 
problems of testing and classification of integrable nonevolutionary equations of orders 4,6 and 10. We believe 
that we have found a few new integrable cases. In all cases we have proven the integrability by presenting 
a bi-Hamiltonian structure, a corresponding recursion operator, a Lax representation or a suitable linearising 
substitution. Our approach can be developed even further, namely a global, i.e. in all orders, classification 
of integrable equations can be achieved (it would require some elements of the Number Theory, such as the 
application of the Lech-Mahler Theorem and it is beyond of the scope of this paper). 
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